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ly-^ ■ We compute the one-loop expressions of two running coupling constants of the Luttinger model. 

f***^ ' The obtained expressions have a nontrivial momentum dependence with Landau poles. The reason 

for the discrepancy between our results and those of other studies, which find that the scaling laws 
£\J . are trivial, is explained. 
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INTRODUCTION 



The Luttinger liquid, namely the one-dimensional system of interacting massless fermions with independent con- 
servation laws for the left- and right-moving species, has been and still is the subject of many theoretical studies in 
^ , Condensed Matter Physics (see, e.g., 0,0,111). This is due to the fact that it has several characteristic properties which 
are very different from those of a Fermi liquid. Indeed, its momentum distribution is continuous at the Fermi energy 
even at zero temperature. Moreover, there are no fermionic quasi-particles in the system because its single-particle 
density of states vanishes as a power law when one approaches the Fermi energy. In addition, the Luttinger liquid 
has collective bosonic charge- and spin-density modes. The equations of motion of the system are solved with the 
help of the method of bosonization jj| . This method allows to replace the complicated fermionic dynamics by another 
one which involves noninteracting bosonic degrees of freedom only. The beta functions associated with the fermionic 
dynamics are believed to be equal to zero 0, S , a result which is commonly regarded as reflecting the triviality of 
the bosonic dynamics. 

The dynamics of the Luttinger liquid may actually be less trivial than expected. It may indeed well be that the 
particles in the system interact also by means of long range forces, as suggested by the following line of argument. It 
is known that the one-dimensional propagator of a free massless fermion has a large anomalous dimension [g, @ ■ This 
large anomalous dimension modifies the singularity structure in momentum space in such a way that the scattering 
amplitude for isolated particles vanishes together with the residue of the propagator on the mass shell 0, 0] . One 
t-H , may therefore say that a given particle moves slower or faster because of the nature, repulsive or attractive, of the 
C ■ interaction rather than because of the relative remoteness or closeness of other particles in the system. Such a 
behaviour clearly suggests that the interactions between massless fermions are able to act over large distances in one 
dimension. The scale parameter controlling this behaviour can only originate from a well-defined cutoff because the 
Luttinger model is known to be classically scale invariant [lCj . 

This line of argument, however, doesn't provide a proof for the existence of long range forces in the Luttinger 
liquid. Our aim in the present paper is to establish in a more reliable way whether such forces exist in the system 
I ■ or not. This is done with the help of a renormalization group study of two effective coupling constants of the model. 
The computations are done in the regularized version of the model with a sharp cutoff. We find that the one-loop 
expressions of the running coupling constants have Landau poles which line up in a scale invariant manner around the 
[ Fermi points and at definite momentum scales resulting from the finiteness of the cutoff. The fact that the momentum 
dependence of these quantities is singular around the Fermi points allows us to conclude that long range forces do 
exist in the Luttinger liquid and so that the dynamics of this system is strongly coupled. 

X 

VERTEX FUNCTIONS 

The regularized Luttinger model with a sharp cutoff of value A is defined by the action S-R eg (A) = So + Si (A), 
where 



S = yjC(k)fc-V-(k)+V4( k ) fc +<M k )] (!) 
describes the free propagation of massless fermions and 

Si(A) = -g [ V;(k + q)V+(k)C(k'-q)V-(k')/'+(A,fe + g)p + (A,fe)p_(A,fe'-g)p_(A,A ; ') (2) 

A,k',q 
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their interactions. This action is the simplest one of all those considered in the context of the so-called g-ology models 
[Tlj . The left- and right-moving spinless particles are represented by the fields ipl) and t/*+) respectively. 
We use the conventions J k = J dk dk / \2ir) 2 and k = (k ,k), where fc and k are the frequency and momentum 
variables respectively. The momenta k± in Eq. are defined by the relation k± — ±kvp — fj, — ik , where vp is 
the Fermi velocity and p = kpvp the Fermi energy (kp is the Fermi momentum). The cutoff functions p± in Eq. 
are Heaviside functions 6; more precisely, p±(A, k) = 0(A — \k =p kp\). We assume in this paper that the absolute 
value of the bare coupling constant g is small. It is important to bear in mind the fact that the expression of So, 
since it contains the linearized form of the dispersion relation, is valid for small values of A only; the limit A — > oo 
is therefore a purely formal limit in the context of the Luttinger model. We use here the method of the equation 
of motion to compute the expressions of the vertex functions needed in the renormalization group study. Although 
these expressions could be obtained by a straightforward application of the Feynman rules, we prefer to follow the 
more compact computational scheme offered by the chosen method because it shows in a more transparent manner 
the effect of the regularization on the dynamics of the model. 

Let us first recall that the generating functional W[rj±,ri±} for the connected 2n-point Green functions G^ 2n ^ 
(n = 1, 2, ...) of the model is defined by 



D[r.]D[^-]DW + ]D[^+]e 



-S R e B (A)-S Sr 



(3) 



with 



S Src = I ^K(k)^ Q (k)+C(k)r/ Q (k)]. 

Jk a=± 



(4) 



One has then (a,; = ±, i = 1, ...,n) 



Gi 2 " ) . Q „ (pi, • • • , p„, pi, • • • , P^) = (2tt) 4 



8 2n W 



5if ai (Pi) • • ' $Vt n (Pn)H*i (Pi) • • • 5v* n (K) 



(5) 



Each Green function is regularized by the product of the cutoff functions corresponding to the external and internal 
lines. However, since all external momenta belong necessarily to the intervals [ikp — A, ±kp + A], the cutoff functions 
corresponding to the external lines can be safely ignored. Let us also recall that the Legendre transform T[i/j±, ip±\ of 
W[?7±) T}±] is the generating functional for the 2n-point vertex functions r^ 2 ™- 1 (n = 1, 2, ...) of the model. Indeed, one 
has [oti = ±, i = 1, n) 



ri 2 ". ) . Qn (pi, ■ • • ,Pn,Pi, • • • , Pn) = (2tt) 4 



S 2n T 



<%! (Pi ) ' ' ' (Pn)^ ai (Pi) • ' ' Stp an (p'„) ' 



(0) 



(4) 



The two effective coupling constants considered in this paper are those related to the four-point vertex functions T 

and rvi., whose expressions are computed below in the one-loop approximation. 

Infinitesimal translations ip a — > tp a + e {a — ±) of the integration variables ip a in Eq. J2J lead to the two regularized 
equations of motion of the model, one for each species of particles. These two equations may be combined into a 
single one, namely, 







(2nfp a 



vUp) 



+g{2nfp a {A lP ) 



k.q 



P-(A,p+ q)p+(A, k - q)p+(A, k) 



fia,+ \ P+ (A, p + q)p- (A, k - q)p- (A, k) 

6 3 



6 3 



8rj- (p + q)6r) + (k - q)^ (k 



(5?7 + (p + q)<5?7_(k - q)^* (k)_ 

w . (7) 



If we let the functional derivative 5/5rj*L(p') act on the regularized equation of motion corresponding to the case 
a = +, we obtain the following regularized Dyson-Schwinger equation: 



= P+ gV 2) (p',p) + <S(p-p')+.9 f 

Jk,q 



G^(p',k,p + q,k q) 



-(2 7 r) 2 G (2) (k,k-q)GV ZJ ( P ',P + q) 



p+ (A, p + q)p- (A, k - q)p- (A, k). 



(8) 
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The action of well-chosen third-order functional derivatives on the regularized equation of motion corresponding to 
the case a = +, Eq. (JTJ, leads to the identities needed to compute the one-loop expressions of the vertex functions 
T^l_ and If we let the functional derivative 5 3 /5r]^(jp')5ri- (r)<5?7* (r') act on this equation of motion and then 

use the regularized Dyson-Schwinger equation, Eq. JSJl, we obtain the following identity: 



= p + GL 4 i(r',p',r,p)- 5 [ 

Jk,q 



G™_ + (r>, k, p', r, k - q, p + q) + (2tt) 2 Gf (p', p + q)GL 4 l (r', k, r, k q 



GL 2) (r', k - q)G^ (p', k, p + q, r) + G^> (k, k - q)G^| (r', p', r, p + q) - G^> (k, r)G% (r', p', k q, p + q) 



,(4) 



-,(2 



( 4 ) W ~' 



.(2), 



(2tt) 4 G (2) (k, r)G^ (r', k - q^ (p', P + q 



(2), 



P+ (A, p + q)p^ (A, k - q)p_ (A, k) . 



(9) 



The one-loop expression of v fl_ is obtained from this identity in the following way. First, we express each connected 
four- or six-point Green function in Eq. in terms of its associated vertex function. Next, we replace each vertex 
function ra lQ2 (p, r, p', r') (oil ^ (X2) that has appeared in the integrand of the equation by its expression at leading 
order in <?, namely, 



r£) Q2 (p,r,p',r') = gS{p + r - p' - r')(<WW + <W-<W) + O 1 



(10) 



(the vertex functions Y__ and r^^ + need not be considered because their expansion in powers of the bare coupling 
constant starts with a power of g larger than one). In a last step, we simplify the integrand by using the expression 
of the regularized Dyson-Schwinger equation at leading order in g, namely, 



P+ + g I G (2) (k)p_(A,fc) + 0( ff 2 ) 



We find that 

r ( f 4 l( P ,r,p') = .g + .g 2 



p+(A,p' + q) p_(A,r + g) P-(A,p + r - p' - q) p+{A,p' + q) 
(p' + Ct)+ (r + q)- {p + r-p' -g)_ (p 1 + q) + 



0(g 3 



(11) 



(12) 



where ryl(p, r, p') is introduced by means of the relation r^_(p, r, p', r') = <5(p + r p' r')r^_ (p, r, p'). 

If we let the functional derivative 6 3 / 5r)*L(r')Sr)+(r)Sr)*i(p') act on the regularized equation of motion corresponding 
to the case a = + , Eq. J7J, and then use the regularized Dyson-Schwinger equation, Eq. (J8J, we obtain another 
identity which is the following: 



= P+ Gf + (r',pi,r,p)-g f 

Jk,q 



GV 6 i+(r', k, p', r, k - q, p + q) + (2tt) 2 Gf (p', p + q)G^ (r', k, r, k - q) (13) 



p+ (A, p + q)p_ (A, k-q)p- (A, k). 



-Gf (r', p + q)G^(p', k, r, k q) + G (T > (k, k - q)G^(r', p', r, p + q) 



The one-loop expression of I++ is obtained from this identity by following the same procedure as for r+! . We find 
that 



lti(p,r,p')=. 9 2 / 



p- (A, q) P-(A, q + r - p') p_ (A, q) P-(A, q + p-p 1 ) 



q- (q + r-p')- 



q- {q + p-p')- 



+ O (g 3 ) , 



(14) 



where r^_(p, r, p') is introduced by means of the relation r^_(p, r, p', r') = <5(p + r — p' — r')T++(p, r, p') 



(4) 



.(4) 



RUNNING COUPLING CONSTANTS 



(4) 

For later convenience, the variables of the vertex functions (cr = ±) are parametrized as 

k = (-iv F {k + A k ),k F + k), 1 = (-iv F {l + A t ), a(k F + I)), 
k' =(-iv F (k' + A k ,),k F + k'), V = (-iv F (l' + A l ,),a{k F +l')), 



(15) 
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where the momentum variables A k , A;, Ay, and A;/ measure the distance from the mass shell. Conservation of 
energy imposes the condition k + A& + I + Ai = k' + Ay + I' + A;< and conservation of momentum the condition 
k+<rl=k'+al' . The combination of these two conditions leads to the additional constraint A^ + A; = Ay + A;/ in 
the case of the vertex function r^ 4 }. The mass shell expression of T^l_ depends on two momenta only (say, k and I) 
because the combination of the two conservation conditions with the condition A^ = A/ = Ay = A/< = gives the 
additional constraints k — k' and I = I' in the case of this vertex function. Performing the integrations in Eqs. i|12|) 
and H14f> . we find that 



rfl(k,i,k') = g. 



/ + _(k,l,k') + 0( 5 3 ), 



ri 4 i(k,i,k') 



Attvf 



/ ++ (k,l,k') + 0(g 3 ), 



(16) 



where 



/+-(k,l,k') 
/++(k,l,k') 



In 



(k' + l) 2 - (k' -1 + A k > - A;) 2 
(k-l) 2 - (k + l + Ak + Atf 
k-k' k' 



hi 



(2A- \k-l\f - {k + l + A k + A[) 2 



(2A - \k' + l\) 2 ~(k' -1 + Ay - A ; ) 2 



k-k' 



(A k 



Ay] 



k'-l 



(A k 



(17) 



The renormalization group method allows us, starting from the truncated expressions of the vertex functions, Eqs. 
(|16fl . I|17l) . to compute the expressions of the running coupling constants we are interested in. This method is usually 
implemented in three steps which, in the present context, are the following. 

(i) In a first step, we introduce two effective coupling constants g^ and g++ by choosing a subtraction point 

(k, 1, k') in the six-dimensional energy-momentum space and setting 



g+- 
g++ 



ri 4 l(k,U'). 



(18) 



The subtraction point is chosen in such a manner that the absolute value of each effective coupling constant is small. 
It has to be noticed that if g^ is already present in the bare action, g ++ is fully generated by the dynamics. 

(ii) Next, we choose a scheme, a trajectory in the energy-momentum space which can be parametrized by a variable, 
denoted by A in the sequel. The trajectory crosses the subtraction point for A = A. The trajectory usually chosen 
is the one associated with an homogeneous rescaling of the variables k, A k , etc., in which case A controls the scale 
dependence. 

(iii) The values of the effective coupling constants g +a (X) (a — ±) corresponding to an arbitrary point of the 
trajectory, that is, to a given value of A, are then computed by integrating the beta functions along this trajectory 
from the subtraction point up to the considered point. The beta functions /3 CT (A) give the rate of variation of the 
effective coupling constants as A increases, namely, 



/?a (A) = A 



dg +a W _ g 2 +a W x df+AX) 



0(g\ 



d\ Attvf dX 

The integration of these quantities leads to the following expression of the running coupling constants: 



g+<r(\) 



g+cr 



l + ^[/ +CT (A)-/W(A)]' 



(19) 



(20) 



with gj ra =g +C j(X). This expression, which gives at once the values of the effective coupling constants at any point of 
the trajectory, may be viewed as resulting from an optimized partial resummation of the full perturbation series of 
the vertex functions. 

Notice that the particular trajectory which defines the scheme does not appear in the expression of Eq. I|2UI) . This 
circumstance allows us to generalize this expression to an arbitrary scheme, namely, 



3+CT (k,l,k') = 



9-1 



l + fe[/+<x(k,l ) k')-/+ < ,(k,I,k')] 



(21) 
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This expression offers the advantage that the effective coupling constants depend only on the coordinates (k, l,k') 
of the point at which they are computed. Its derivation is based upon the fact that the rate of change of the 
renormalized coupling constants caused by an infinitesimal step in the energy-momentum space has a component 
along each direction of this space. The expression of such a component is 



/V(A) 



dq 



Attvf 



(q) q df +CT (q) 



dq 



O 



(9lo 



(22) 



where the symbol q represents a component of (k, 1, k'). The integration of these "directional" beta functions along the 
chosen trajectory leads to the expression of Eq. 121|) which, interestingly enough, may also be viewed as originating 
from an optimized partial resummation of the Bethe-Salpeter equation. 

A mass shell expression of g^ which is valid in the particle-particle (p-p) channel is obtained by using Eq. 12Ufl in 

the case of the linear trajectory fc(A) = /(A) = AA (|A| < 1). We find that 



9 P +- P W 



1 + 



'(A) 



hi 



(SSI)-* (SGI) 



(23) 



This equation clearly diverges and the Landau pole is found at 



Ar p = 



1 + A - (1 - A)e 



(24) 



1 + A+(1-A)e s +- 



The effective interaction between the particles has obviously not the same sign in the particle-hole channel as in the 

particle-particle channel. It is easy to verify that the function /_| , Eq. l|17Jl . the one used to define the effective 

coupling constant , Eq. (|18fl . changes sign if the momentum I is tranformed into —I; a mass shell expression of 

which is valid in the particle-hole (p-h) channel can therefore be obtained by using Eq. (|20|l in the case of the 

linear trajectory fc(A) = —/(A) = AA (|A| < 1). The result is 



9 P + --(X) = 



9 P +-W 



1 - 



In 



i-|A| 



(25) 



This equation diverges too and the Landau pole is at 



1 + A-(1-A)e 9 +- 



(26) 



1 + A + (1-A)e 9 +- 
We use now Eq. I|21|) to obtain g^ on the mass shell, 



9+-(k,l) 



1 



9+- 



a+- 



In 



(2A-|fc-i|) 2 -(fc+Q 2 



(2A-|fe+i|) 2 -(fe-0 2 



-In 



(2A-|fc-i|) 2 -(fc+Q 2 



(2A-|fc+;|) 2 -(fc-z) 2 



(27) 



Notice that the first term on the right-hand side of the first equation in Eq. H17|) is vanishing on the mass shell. As a 
result, the scale dependence, i.e. the deviation from a scale invariant result, arises from the presence of the cutoff A 
only. The Landau poles corresponding to the divergences of this expression belong to the line of the plane (k, I), 



Attvf 
9+- 



ln 



(2A- \k-l\) 2 - {k + lf 
(2A- \k + l\) 2 - {k-iy 



In 



{2K~\k-l\f - {k + lf 

(2A- |fc + ij) 2 - {k-Ty 



(28) 



This line stays away from the Fermi points for any generic choice of the bare coupling constant, that is, in the absence 

of fine tuning. The behaviour of g^ that is valid slightly away from the mass shell is similar to the one of g^ on 

the mass shell and the corresponding equations can be obtained from Eqs. J23|)-(|2E1) by adding terms of the order 
O (A) on the right-hand sides. 
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The coupling strength g ++ is always vanishing at the order of one loop on mass shell. Using Eq. 121fl again, we 
obtain 

(29) 



k-k'+HA k -A k ,) fe'-i+i(A k ,-A,) k-k' + ±(A k -A k ,) + k>-l+±(A k ,-A,) 



4-7TVF 

off-shell. The Landau poles are now spread over a five-dimensional hypersurface whose implicit equation is 
A-kvf k-k 1 k'-l k-k' k' — l 



~9++ fc-fc' + i(A fe -A fe ,) fc'-Z + ±(A fe , -A,) fc-fc' + i(A fe -A fc k'-l + U& k 



(30) 



This hypersurface approaches the Fermi point of the right-moving particles in the limit Afe, A;, and Afe' — > for any 
value of g++. 

Several comments are in order at this point. We start with the scale dependence of the coupling constants and 

g++- The latter is scale invariant according to Eq. 1(29(1 and its beta function is vanishing in the conventional schemes 

where the variables k, A/., etc. are scaled homogeneously. The scale invariance of g-\ is broken by the cutoff, as 

shown by Eq. 1)27(1 . The limit A — > oo carried out in a given finite order of the perturbation expansion suppresses this 
breakdown and one would naively expect that the renormalized theory recovers the classical scale invariance. But 
the partial infinite resummation of the perturbation expansion, carried out by the generalized renormalization group 
scheme, indicates that the interactions become strong within certain kinematical regions because of the accumulation 
of large contributions of higher orders in the perturbation expansion. These regions are characterized by the scale of 
the cutoff. This prevents us from removing the cutoff and restoring the naive, classical scale invariance at the same 
time. 

Another remark is about the nature of the singularities of the effective coupling strengths. The divergence in Eqs. 
(|27|l and l(29() does not necesserily mean that the theory becomes strongly coupled. In fact, the effective interaction 
strengths should be defined by scattering amplitudes or other physical observables which reflect the interactions among 
elementary or dressed excitations. Such a construction always involves a finite resolution in the energy-momentum 
space which can be imagined as coming from a smearing of some Green functions in the kinematical space. As 

long as the singularities are integrable, as is the case of <?_| in Eq. (|2*7jl . such a smeared coupling strength remains 

finite despite the diverging numerical values the effective coupling constant may take in the vicinity of the Landau 
poles. But the non-integrable singularity of g ++ appearing in Eq. (|29|l produces a diverging coupling strength for an 
arbitrary localisation of the constituents in the energy-momentum space. 

We close this Section by pointing out the way the apparent contradiction between our results and the well-known 
vanishing of the usual beta function can be resolved. 

(i) The beta function computed in Refs. 0, @ is based upon a particular choice of the subtraction point, namely 
Po = r = p'o = ?o = (A fe = -k, A ; = -I, A k > = -k', AJ = -I') and k = -I = —k' = -l'/3. The special feature of 

this point is that the radiative corrections /_| and /++, given by Eqs. (|17|l . are vanishing. This subtraction point 

together with the usual, homogeneous scaling implies vanishing beta functions. But the generalized scheme used 
in this work picks up the radiative corrections arising from an arbitrary displacement in the kinematical space and 
reveals the non-naturalness of this cancellation. 

(ii) Another argument, the one used in Ref. 0, rests on the study of the behaviour of the bare coupling constant 
in the blocking procedure. The bare coupling constant was defined in that work at the Fermi point and the possible 
energy-momentum dependence was treated perturbatively. The point is that the leading order result, the O (<? 2 ) 
coupling constant at the Fermi point, agrees with Eq. ((12(1 . The singular energy-momentum dependence generated 
by the partial resummation in Eqs. l(2"7| and (|2"9"|) suggests that the dependence of the coupling constant on the 
energy-momentum can not be taken into account by expanding around the Fermi point and the simplicity of the 
scenario spelled out in Ref. appears misleading. 

(iii) A more formal argument has also been used to justify the triviality of the scaling laws. This argument is based 
upon the (supposed) equivalence between the behaviour of the running coupling constant in two different limits, the 
one where the cutoff is removed and the one where each energy-momentum variable k tends towards zero. Indeed, 
since these variables and A are the only scales (fcp- is absent from the beta functions), the dimensionless running 
coupling constants should not depend on them separately but rather on the ratios k/A. Consequently, the limit 
k — > is supposed to coincide with the limit A — > oo where the scale dependence drops, as one can see by inspecting 
Eqs . ((17(1 . The fallacy of this argument is usually shown by recalling the phenomenon of dimension transmutation 
[12j. This phenomenon is found in asymptotically free, classically scale invariant models where the cutoff generates a 
scale which remains finite even in a renormalized theory where the cutoff is already removed. The Luttinger model 
offers a simpler mechanism to break the scale invariance by its Landau poles, as pointed out after Eq. I(3U|) above. 
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SUMMARY AND CONCLUSION 



Our aim in this paper was to establish the existence of long range forces in the Luttinger liquid in a reliable manner. 
This has been done with the help of a renormalization group study of two effective coupling constants of the model. 
The computations have been done in the one-loop order. We have shown that the running coupling constants have 
Landau poles in the vicinity of the Fermi points for arbitrarily weak bare interactions. The singular momentum 
dependence of the effective coupling constants in the vicinity of the Fermi points shows that the interactions between 
the fermionic excitations are of long range. 

The problem of the breakdown of the perturbation expansion in the vicinity of the Fermi points is difficult to solve. 
It is known that the correspondence between the fermionic and bosonic versions of the model can be established at any 
finite order of the perturbation expansion in powers of the particle mass or interaction strength. However, since the 
strong interactions close to the Fermi points modify the single particle dispersion relation in a fundamental manner, 
the perturbative approach may not be sufficient to prove this correspondence. We believe that a more systematic 
study of the scaling laws established in both versions of the theory is needed to clarify this question. 

The partial resummation, performed by the integration of the beta functions, remains consistent as long as the 
effective coupling strengths are weak. When the one- loop beta functions predict large values for the coupling constants, 
the computation is no longer reliable and one usually assumes that the true dynamics is indeed strongly coupled in 
terms of the degrees of freedom used. A similar problem, namely the one of the soft gluons, blocks the way towards 
the analytical understanding of Hadronic Physics. We believe that any progress towards more definite conclusions 
about the dynamics of the Luttinger model requires some non-perturbative method in order to treat the strong, long 
range interactions, just as in the case of QCD. 

We wish to acknowledge useful discussions with D. Cabra and J. Solyom and to thank K. Kikoin for bringing Ref. 
to our attention as well as the Referee for helping us to improve the content of the paper. 
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